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We investigated the influence of efficacy of synaptic interaction on firing synchronization in exci- 
tatory neuronal networks. We found spike death phenomena, namely, the state of neurons transits 
from limit cycle to fixed point or transient state. The phenomena occur under the perturbation 
of excitatory synaptic interaction that has a high efficacy. We showed that the decrease of synap- 
tic current results in spike death through depressing the feedback of sodium ionic current. In the 
networks with spike death property the degree of synchronization is lower and unsensitive to the 
heterogeneity of neurons. The mechanism of the influence is that the transition of neuron state 
disrupts the adjustment of the rhythm of neuron oscillation and prevents further increase of firing 
synchronization. 
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I. INTRODUCTION 

Synchronization of neural activity appears in different 
parts of the mammahan cerebral cortex [l| , and underhes 
different neural processes in both normal and anomalous 
brain functions Q. It has been suggested that synchro- 
nization plays a vital role in information processing in the 
brain, e.g., processing information from different sensory 
systems to form a coherent and unified perception of the 
external world [l|, y, y, Ij, S S]- On the other hand, syn- 
chronization has been detected in pathological conditions 
such as Parkinson's disease 0, @] . And epileptic seizures 
have long been considered resulting from excessive syn- 
chronized brain activity [9|, though some recent studies 
suggest that this picture may be an over-simplification 
[lOl . Ill| . Therefore understanding the mechanisms of syn- 
chronization may be a critical step in elucidating how 
neural systems work [il[. It has stimulated a great deal 
of theoretical and numerical works, such as the studies 
on the effects of the topological properties of underlying 
networks [ij, llj, H-llal a nd the dynamical properties of 



synaptic coupling [T^, [Tjl 

It was recently shown that the response time of synap- 
tic couplings influences the stability of synchronized os- 
cillation in the nonlocally coupled Hodgkin-Huxley (HH) 
equations [l6| . If the response time of synaptic coupling 
is slower, synchronized activity of the systems is instable 
for excitatory coupling. However, the underlying dynam- 
ical mechanism of the influence is not clear. In experi- 
mental studies [l^ , it has been suggested that the gener- 
ation of prolonged epileptiform neuronal synchronization 
is favored by lower efficacy of synaptic transmission. The 
numerical studies [19[ in a detailed computational model 
revealed that seizure-like activity occurs when the exci- 
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tatory synapses are weakened, and the results were con- 
firmed experimentally in mouse neocortical slices. Ac- 
cording to the common accepted assumption that syn- 
chronization of neuronal activity underlies seizures, the 
dynamical mechanism of synchronization may be useful 
for understanding the way the biological neural system 
works. 

In this work, we numerically investigated the dynam- 
ical mechanism underlying the influence of synaptic ef- 
ficacy on firing synchronization in HH neuron networks. 
To do this, we first studied the dynamics of the response 
of HH neuron to excitatory synaptic current. When the 
efficacy of synapse is low, namely, strength is weak and 
duration is short, the limit cycle is stable to the pertur- 
bation of the synaptic current. When synaptic efficacy 
is high, synaptic current can induce the transition of the 
neurons from limit cycle to fixed point or transient state. 
The transition is determined by dynamics of neuron's 
ionic channel. The decrease of synaptic current depresses 
the feedback of sodium ionic current which is responsible 
for the initiation of the spike. For simplicity we will refer 
to the transitions as spike death. 

In neuronal networks, synaptic input of a neuron is the 
accumulation of the currents received from all presynap- 
tic neurons. When the coherence of firing time of neurons 
is enhanced by the excitatory interaction, the synaptic 
input of neurons transforms from the fluctuating wave- 
form into the pulse shape like the signal produced by 
one synapse. If synaptic efficacy is high, the input sig- 
nal can induce spike death of the neuron. Then spike 
death disorders the adjustment of the rhythm of neurons 
and prevents neurons from firing spikes synchronously. 
In contrast, for synapses of lower efficacy, the duration 
of synaptic current is too short to induce spike death 
of neurons. Additionally, the firing synchronization is 
different from synchronous activity of oscillators for the 
existence of the transitions of neuron's state. 

The paper is organized as follows. The HH neuron 



model and the synaptic coupling were introduced in Sec. 
II. The response of a HH neuron to synaptic current was 
investigated in Sec. III. The influence of the dynamics 
of neurons on firing synchrony was shown in Sec. IV. 
Discussion and conclusion were given in Sec. V. 



II. MODEL 

To investigate the dynamics of a neuron under the per- 
turbation of synaptic stimulus, we adopted a system con- 
sisting of a HH neuron and a synapse. The HH neuron 
was originally proposed for the giant axon in a squid 
[20[ . It serves as a paradigm for the spiking neuron mod- 
els that based on nonlinear conductances of ion channels. 
The model describes the evolution of the membrane po- 
tential V{t) and can be written as 



fixed point becomes unstable. The values of bifurcation 
points are /g ~ 6.2 //A/cm^, /i w 9.8 /iA/cm^ and I2 ~ 
154^A/cm2 m. 

We adopted the synaptic current I^yn described by an 
alpha function [22| . The alpha function synapse is a 
phenomenological model based on an approximate cor- 
respondence of the time course of the waveform to physi- 
ological recordings of the postsynaptic response [2J| . The 
equation of the synapse is like: 



Isynit) = -gsyna(t - tin)(y{t) - Egyn), (10) 



with 



a{t) = (i/T)exp(-t/T)e(t). 



(11) 



ii ^~ -^ion I ^stini ~r J-syn^ 



(1) 



where Lnn is the ionic current, /, 



is the external cur- 



Where r is the characteristic time of the interaction, 8(t) 
is the Heaviside step function, and ij;„ is the beginning 
time of the synaptic interaction, i.e., the firing time of 
the presynaptic neuron (all delays are neglected). The 
synaptic effect is traditionally classified as excitatory or 
inhibitory depending on the value of Esyn- Here, we 



rent and Isyn is the synaptic current. The ionic current took Esyn = 30 mV for excitatory synapses, and —80 
describes the ion channel on the membrane and is defined mV for inhibitory ones. Eq. pT|) yields pulses with 



the maximum value of 



aL T — T-j 



T and with the 



Iton = -gNam^HV - E^a) - dK-n^iV - Ek) - gi{V - El), 

(2) 
where gNa, gK and gi are the maximum conductances for 
the sodium, potassium and leak currents, -Eato, Ek and 
El are the corresponding reversal potentials, m and h are 
the activation and inactivation variables of the sodium 
current and n is the activation variable of the potassium 
current. The gating variables y = m, h, n satisfy the 
differential equation 



dyjt) 

dt 



^ay[l-y{t)]-Pyy{t), 



(3) 



with nonlinear functions ay and (3y given by 



a™ - 0.1(F + 40)/{l-exp[-(T/ + 40)/10]}, (4) 

(3m = 4exp[-(F + 65)/18], (5) 

ah = 0.07cxp[-(y + 65)/20], (6) 

(3h = l/{l + exp[-(T/ + 35)/10]}, (7) 

a„ = 0.01(T/ + 55)/{l-exp[-(F-F55)/10]}, (8) 

/3„ = 0.125cxph(F + 65)/80]. (9) 

The parameter values are E^a ~ 50 mV, Ek = — 77 mV, 
El = -54.4 mV, gNa = 120 mS/cm^, off = 36 mS/cm^, 
gi = 0.3 mS/cm^, and C = 1 fiF/cm'^ MH^. 

The external current IsUm determines the firing rate 
of the neuron. In the absence of synaptic coupling Igym 
the HH neuron has the following bifurcation diagram as 
a function of IgUm ■ In the parameter regions IsUm < ^o 
and IsUm > I2 fixed point is the global attractor. For 
-^0 < Istirn < Ii the neuron possesses coexisting stable 
attractors, the fixed point and the limit cycle which are 
separated by an unstable limit cycle. For /i < IsUm < I2 



half- width of 2.45r [25[. So r characterizes the dura- 
tion of the synaptic interaction. For the alpha function 
synapse, we equated the synaptic efficacy to the maxi- 
mum synaptic conductances gsyn and the characteristic 
time T. The high efficacy means that synaptic current 
possesses strong strength and long duration. 



III. SPIKE DEATH OF NEURON 

We focused on the dynamics of the system in the pa- 
rameter region near the bifurcation point /i . Firstly, we 
studied the response of bistable neuron (/q < IsUm < h ) 
to excitatory synaptic current. In simulations, firing was 
identified as membrane potential V is over 20 mV. When 
the neuron fired a spike, we triggered a pulse of synaptic 
current into it. We observed two types of dynamics of the 
response, which depended on the efficacy of synapse. For 
slow response time and strong synaptic strength, the neu- 
ron transited from the limit cycle to the fixed point. The 
transition of neuron's state is shown in Figs. [1] In Fig. 
[IJa) the response of the neuron to the synaptic current 
is represented by the membrane potential V. In the fig- 
ure, the synaptic current is illustrated by the pulse added 
on the external current. One can see that the periodic 
firing eliminated after the synaptic current was injected. 
The value of membrane potential tended to -61.15 mV 
through subthreshold oscillation. In Fig. [TJb) the tran- 
sition of neuron's state is shown in a three-dimensional 
space {V, h, m) which is a projection of the phase space 
{V , h, 771, n). The trajectory left the limit cycle and 
was attracted into the basin of the fixed point. At last, 
through transient process, the trajectory stopped at the 
fixed point which is indicated by the dashed lines in the 
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FIG. 1: (Color online) (a) The response of bistable HH neuron 
to excitatory synaptic current. The symbol I denotes the 
sum of the external current Isum and the synaptic current 
Isyn- (b) The corresponding phase portrait. The values of 



parameters are r 
/iA/cm^ . 



2 msec, gsyn = 1 mS/cm and Is, 



figure. The fixed point was {V,h,m,n)={-60.15, 0.423, 
0.092, 0.394) as the parameter IsUm = 8.5 /lA/cm^. On 
the other hand, for systems with quick response time and 
weak synaptic strength, the transition did not occur and 
the trajectory was attracted back the limit cycle from 
weak perturbation. 

When the external signal IsUm is larger than but near 
/i, the neuron possesses a stable limit cycle and an un- 
stable fixed point. In this case there were also two types 
of dynamics of response. For high synaptic efficacy, the 
neuron exhibited a transient behavior when it received 
the synaptic current. Fig. [2] (a) shows that the mem- 
brane potential responded to the synaptic current with a 
transient subthreshold oscillation. The subthreshold os- 
cillation interrupted the periodic firing of the neuron. In 
phase space the transient behavior was a motion around 
the unstable fixed point. This is explicitly shown in the 
corresponding 3D phase space {V, h, m). In Fig. [D^b) 
one can see that the trajectory of the neuron left the limit 
cycle and transiently moved around the unstable fixed 



FIG. 2: (Color online) (a) The response of the neuron with 
IsUm = 12.5 /iA/cm^ to excitatory synaptic current, (b) The 
corresponding phase portrait. The synaptic parameters are 
the same as Fig. [1] 



point which is indicated by the dashed lines. The unsta- 
ble fixed point was {V,h,m,n)={-58.704, 0.374, 0.108, 
0.417) as the parameter IsUm = 12.5 ^A/cm^. On the 
other hand, for low synaptic efficacy, the coupling can- 
not induce the transient motion around the unstable fixed 
point. Like bistable neurons, the trajectory returned to 
limit cycle from weak perturbation. 

Fig. [3] shows the boundary between the two types 
of dynamics on the parameter plane gsyn vs r. For 
the bistable neuron of the external current IsUm = 8.5 
/lA/cm^, the boundary between the two types of dynam- 
ics is represented by squares. Above the curve, neurons 
responded to synaptic currents with transitions between 
attractors. Inversely, the limit cycle of neurons was sta- 
ble to the perturbation. For the neuron of the external 
current IgUme = 12.5 /iA/cm^, the boundary is shown by 
circles. It is notable that stronger strength and longer 
duration of synaptic current were needed by the transi- 
tions of neuron's state. Based on numerical simulations, 
we obtained that the transitions of neuron's state can 
be induced by synaptic current in the parameter region 
6.2 ^A/cm^ < Istim < 28.8 ^A/cm^. Above the upper 
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FIG. 3: (Color online) The dynamics diagram of minimal 
synaptic intensity gsi/n vs the characteristic time r for the neu- 
rons with the external current IsHm = 8.5 ^A/cm^ (squares) 
or Istim ~ 12.5 /iA/cm^ (circles). 



FIG. 4: (Color online) The change of the activity of HH neu- 
ron as external current decreases. The external current is 
reduced from 20 /xA/cm to 8.5 /iA/cm with the slop (a) -1 
/iA/(cm^ msec), (b) -2 /iA/(cm'^ msec), or (c) -0.5 /^A/(cm^ 
msec) . 



boundary of the region, the transient motion around the 
unstable fixed point cannot be induced by the synaptic 
current. 

In the following, we referred to the phenomena shown 
in both Fig. [l]and Fig. [2] as spike death. Different from 
oscillator death which is the quenching of oscillation of 
coupled systems, spike death is the behavior of a single 
neuron. And spike death includes both the transition 
between stable attractors and the transient behavior. 

Next we gave a qualitative interpretation of the spike 
death phenomena. Although the synaptic current con- 
sists of a rise and a decay stage, we want to show that the 
decrease of the synaptic current induces the transition of 
neuronal activity. We illustrated the role of the decreased 
current in Fig. 01 In the simulations, external currents 
linearly decreased from 20 /iA/cm^ to 8.5 fiA/cm^, and 
synaptic current was absent. The time of beginning to 
decrease and the different slopes of the decrease current 
were chosen to represent the distinct duration and rate 
of decay of the current. In Fig. H] (a), the activity of the 
neuron transited from the periodic firing to silent state 
through transient subthreshold oscillation. The decrease 
of current induced the spike death of the neuron. In Fig. 
m (b), the decay of the current occurred in the refractory 
period of the neuron. The decreased current did not de- 
press the next spike. This is similar to the scenario of 
the synaptic current of short duration. So spike death 
requires that the decrease of current occurs at the end of 
refractory period and the stage of the initiation of next 
spike. In Fig. S] (c), the external current decreased with 
the slope -0.5 ^A/(cm^ msec). Comparing with Fig. |4] 
(a), the rate of decrease was small. Although the current 
decreased at the stage of initiating a spike, it did not de- 
press the firing of the neuron. Thus spike death requires 



that the rate of decrease of current is large at the stage 
of initiating spikes. 

To interpret the effect of the decrease of current on 
the oscillation of neurons, we reviewed the generation 
of spikes. For a fixed value of membrane potential 
V, the variable y (= m, h, n) approaches the value 
yo{y) = C(y{y)/[oiy{V) + Pyi^)] with the time constant 
Tylv) = [ay{V) + l3y{V)]~\ The variable mo(F) in- 
creases with V , and the corresponding time scale t„i is 
littler than t^ and r„. If external current injects into 
the cell and raises the membrane potential V, the con- 
ductance of sodium channels gNa'm-'^h increases due to 
increasing m. Then sodium ions flow into the cell and 
raise membrane potential even further. If this positive 
feedback is large enough, a spike is initiated [2l|. When 
external current decreases at the onset of the generation 
of a spike, the raise of membrane potential slows down. 
Then the increase of m is slowed down. So the positive 
feedback is weakened. If the current decreases quickly, 
the spike of neuron can be depressed. Therefore, for de- 
pressing the positive feedback of sodium ionic current, 
the decrease of current must occur at the onset of gener- 
ation of spike, and the current must decrease in a large 
rate. This interpretation is consistent with above simu- 
lated results that the strong strength and the long dura- 
tion of synaptic current are necessary for spike death. 



IV. INFLUENCE ON FIRING 
SYNCHRONIZATION 

Next we investigated the influence of spike death on 
the firing synchronization in neuronal networks. We con- 
sidered a directed random network consisting of A^ non- 



identical HH neurons. The network was generated as fol- 
lowing: With a probability p, we connected each of the 
probable directed couplings (such as the one from the 
jth neuron to the zth neuron). The network is described 
by adjacency matrix {a*-'} which entry a'-' is equal to 
1 when the coupling from j to i exists, and zero other- 
wise. In the directed network, a*^ can be not equal to 
a-'* , and the signal travels in only the direction from j to 
i if a^^=l and a^^~0. The signal received by the neuron i 
is the accumulation of all input synaptic currents, which 
is defined as 






N 



J = l 



a^^g,yMt-nnWKt)-E%n), (12) 



where (f is the rescaled factor which equals the number 
of inputting synapses of neuron i (i = 1, . . . , N). t^^ is 
the latest firing time of presynaptic neuron j. E^Jy^ is the 
reverse potential of the synapse connecting j to i. 

To study the global behavior of neuronal networks we 
computed the average activity V{t) = (l/N) J2i=i ^*(0 
of the network. The amplitude of average activity can 
intuitionally reveal the coherence of activity of neurons, 
which is defined as [131 
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FIG. 5: (color online) (a) The amplitude of average activity a 
versus the fraction of excitatory neuron fexc- Squares (cycles) 
represented the relations in networks with t—2 msec (1 msec). 
The error bars were the standard deviation across 20 realiza- 
tions, (b) The relation between average cross correlation K 
and the fraction of excitatory neuron fexc- 



a' = 



1 



T2-T1 



[{v{t)), - v{t)rdt, (13) 



Ti 



where the angle brackets denote temporal average over 
the integration interval. In this work wc studied the fir- 
ing synchrony for the reason that neuronal states may 
transit away from the limit cycle and the synchroniza- 
tion of oscillators will be disrupted. Here we focused on 
the coherence of firing time of neurons. We adopted the 
average cross correlation of firing time of neurons |26l . |27[ 
to quantify the degree of firing synchronization. Average 
cross correlation is obtained by averaging the pair coher- 
ence Kij{'y) between neuron i and j, i.e.. 



K 



1 



N 



N 



N{N -1) 






(14) 



The pair coherence Kij('y) is defined as 



K^Jil) 






(15) 



which is measured by the cross correlation of spike trains 
at zero time lag within a time bin 7. To transform the 
neuronal activity into spike train, the interval T2 — Ti is 
divided into k bins of 7 = 1 msec. Then spike trains of 
neurons i and j are given by X{1) = or 1 and Y{1) = 
or 1 (l = 1, . . . , fc), where 1 represents a spike generates 
in the bin and otherwise. 

We numerically investigated the collective activity of 
networks versus the fraction of excitatory neurons fexc 
in the networks. In general, excitatory synapses tend to 



synchronize the activity of neurons in networks. How- 
ever, it was shown that in neural systems different types 
of dynamics of response of neurons to couplings produce 
different synchrony property [23. Here we studied the 
influence of the new type of response, spike death, on the 
formation of synchrony as the excitatory interactions in- 
crease. In simulations, the networks consisted of 1000 
neurons and connecting probability was p^O.Ol. The 
synaptic conductance took the value gsyn=^ mS/cm^. 
As an example, we used the synapses with the character- 
istic time T=2 msec or t=1 msec to generate neuronal 
networks with or without spike death, respectively. To 
ensure that the neurons had nonidentical properties, the 
external currents /^jj^ were in (8.0, 12.0) /iA/cm^ and 
were generated at random. 

In Fig. HJa) we plotted the amplitude of average ac- 
tivity a versus the fraction f^xc of excitatory neurons. 
Squares represents the results obtained in networks with 
T=2 msec, and circles represents the results with t=1 
msec. In Fig. [HJb) we did the same for average cross cor- 
relation K . The relations showed that the coherence of 
activity increases with the fraction of excitatory neurons. 
It is notable that the networks with r=2 msec, i.e., with 
the spike death property, had obviously lower values of 
a and K than networks with t — 1 msec. So the degree 
of synchrony in networks with spike death property was 
obviously lower than networks without the property. 

The difference of the degree of synchrony between two 
types of networks is shown more intuitionally by the aver- 
age activity in Fig. [6l We made neurons in the networks 
beginning to oscillate with a high degree of coherence. 
In simulations the first firing time of neurons randomly 
distributed in (0, 5) msec. When /exc=0.5, the aver- 
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FIG. 6: Average membrane potential of the networks with 
T — 2 msec (left) and t — 1 msec (right). The fraction of 
excitatory neurons is fexc = 0.5, 0.8, 1.0 from top to bottom. 
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FIG. 7: (a) The input synaptic current of a randomly chosen 
neuron from the network with r = 2 msec, f^xc ~ 0.5. (b) 
The input synapse current when fexc increased to 0.8. (c) 
The histogram of spike death events in the network with r = 
2msec and fexc = 0.8. 



age activity of the network with r = 2 msec decreased 
obviously and the network changed quickly to oscillat- 
ing randomly, as shown in Fig. [6l^a). In the network 
without spike death property, coherent oscillation |12l | 
appeared, as shown in Fig. ^h). When the fraction of 
excitatory neurons increased, the excitatory interaction 
enhanced the average activity in the networks both with 
and without the spike death property, while the average 
activity in the latter was obviously larger. Furthermore, 
the frequency of the average activity V in networks with 
T = 1msec was higher than the frequency in networks 
with r=2 msec. 

The synchrony properties can be qualitatively ex- 
plained by spike death. In networks, if the fraction of ex- 



citatory neuron is increased, the coherence of firing time 
of neurons is enhanced. Then the input synaptic current 
of a neuron may transform from fluctuating waveform 
into a smooth pulse. If the pulse of the accumulated 
synaptic current has long enough duration to depress the 
next spike, the synaptic current may lead to spike death. 
Therefore the synaptic current disorders the adjustment 
of the rhythm of neurons' firing and prevents synchro- 
nization. To demonstrate this, we showed the synap- 
tic currents and spike death events in Fig. 7. In Fig. 
7(a) and 7(b) we plotted the input synaptic currents of 
a neuron, which was randomly chosen in networks with 
/e£cc=0.5 and /ea;c=0.8 respectively. The transform from 
fluctuate waveform to smooth pulse was observed. In 
simulations we took the peak of subthreshold oscillation 
of a neuron as a spike death event. Fig. 7(c) shows the 
histogram of frequency of spike death events, i.e., the 
number of spike death events in 1 msec. The spike death 
events periodically appeared with the same rhythm as 
the average activity of the network. In contrast, spike 
death events cannot be obtained in the networks of low 
synaptic efficacy. Thus the spike death can explain the 
above synchrony property. 

If network consists of identical neurons, the influ- 
ence of spike death also exists and is more remarkable. 
Fig. 8(a) shows the relation of K to fexc of the net- 
works in which the external currents of all neurons was 
-fsiim = 10.0/iA/cm^. One can see that the values of K of 
the networks with r=2 msec (squares) was similar to Fig. 
5. For networks with t=1 msec (circles), however, the 
value of K tended to 1 as fexc increases. We calculated 
the degree of synchrony as a function of the heterogene- 
ity of neurons in purely excitatory networks. In simula- 
tions the value of external currents IsUm distributed in 
the region (10.0 - O.Sw, 10.0 + 0.5w) /xA/cm^. Fig. 8(b) 
shows the relation between K and the width w of the 
parameter region of IsUm ■ For networks with t = 2msec 
(squares), the degree of synchrony was unsensitive to the 
heterogeneity of neurons. In contrast, the degree of syn- 
chrony in networks with r = Imscc (circles) remarkably 
increased and tended to 1 as the heterogeneity of neurons 
decreased. For networks with r := 1msec, the relation be- 
tween K and the width w of the parameter region was 
fitted to the first order exponential decay curve. The fit- 
ted curve (dashed line) is K = Aexp{—w/B) + Kq with 
Ko = 0.362 ± 0.005, A = 0.595 ± 0.007 and B = 1.017 
± 0.030. The remarkable difference between two kinds of 
networks shows that spike death can effectively prevent 
the firing synchronization in neuronal networks. 



V. DISCUSSION AND CONCLUSION 

We have used HH neuron networks to investigate the 
dynamical origin of the influence of synaptic efficacy on 
the firing synchronization. A new dynamics of response 
of neurons to coupling, spike death, was suggested as 
a possible mechanism underlying the influence. When 
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FIG. 8: (Color online) (a): The average cross correlation K 
versus the fraction of excitatory neurons for the networks con- 
sisting of identical neurons. Squares (circles) represent the 
relation in networks with r = 2nisec (Inisec). The error bars 
were the standard deviation across 20 realizations, (b): The 
average cross correlation K changes with the width w of pa- 
rameter region of latim ■ Dashed line is the fitted exponential 
curve. 



the firing time of neurons is so coherent that synaptic 
currents have a pulse waveform in excitatory networks, 
synaptic current induces the transition of neuron's state 
from limit cycle to fixed point or transient state. The 
transitions disrupt the adjustment of rhythm of neurons' 
oscillations and prevent further increase of firing synchro- 
nization. 

We studied the dynamics of the HH neuron responding 
to the excitatory synaptic perturbation. We numerically 
demonstrated that the synapse of high efficacy, i.e., large 
characteristic time and strong strength, induces spike 
death of the neuron. For bistable neurons, spike death 
means the transition from limit cycle to fixed point. For 
the neuron with unstable fixed point, spike death means 
the transition from limit cycle to transient state. The 
transient state is the motion around the unstable fixed 
point in phase space. Spike death of neurons results 
from the decrease of synaptic current, which depresses 
the feedback of sodium ionic current at the stage of ini- 
tiating a spike. 

We demonstrated the influence of spike death on the 
degree of firing synchronization, fn simulations wc con- 
sidered the networks with or without the spike death 
property, which were generated using synapses of charac- 
teristic time r=2 msec and r=l msec, respectively. The 
degree of synchrony of the former was lower. This is con- 
sistent with results of [lg| that synchronous state is not 
stable for the excitation coupling of slow response time. 
However, we also showed that for slow response time, the 
degree of synchrony increased with the fraction of excita- 



tory, and oscillation rate of whole network slowed down. 
Our main results are that in the network with t=2 msec, 
the spike death events were observed. And spike death 
can explain the mechanism of preventing the raise of the 
degree of synchrony. Related synchrony proper ties was 
found also in weakly coupled HH neurons [22, [2^ . In the 
case of weak coupling, the phase of neuron was perturbed 
by couplings, but the oscillation of neuron was not de- 
stroyed. In contrast, the new dynamical mechanism we 
suggested is proper for strong coupling and underlies the 
synchrony of the interrupted oscillations. Additionally, 
it is notable that, for the existence of spike death, the 
firing synchronization of neuronal networks is different 
from the usual oscillator synchronization in which each 
oscillator is stable to perturbations [23 . 

Our work relates to that of Drover et al. [30|. In 
the elegant work, using a simplified neuron model of two 
variables, they proposed a mechanism for slowing firing 
down. They found that slow decay synaptic variable in- 
duces that trajectory is attracted toward the unstable 
fixed point of the simplified model. This is similar with 
the transient behavior of HH neuron we proposed here. 
However, with their mechanism, synaptic excitation is 
strongly synchronizing in networks in contrast with that 
spike death prevents synchrony. In neural networks, col- 
lective behaviors sensitively depend on the intrinsic dy- 
namics of neurons [30J, and many types of response of 
neuron to synaptic coupling may exist [23. It is inter- 
esting to make further studies on the relations among 
different responses and their effects on the collective be- 
haviors of networks. 

The variability of strength of synapses was not involved 
in the present investigation. The synaptic strength is 
often affected by the activity of neurons through synaptic 
plasticity [3l| and synaptic adaption [S^]. The effect of 
changes of synaptic strength will be studied elsewhere. 

As mentioned above, the phenomena that low efficacy 
of synapses favors the generation of neuronal synchro- 
nization underlying seizure were obtained in experiments 
and numerical simulations [l8J, [l9J. Here we equated 
the strength and especially the characteristic time of 
synapses to the synaptic efficacy, and studied the mecha- 
nism by which synapses influence firing synchronization. 
The mechanism of the influence may have potentially val- 
ues for understanding the way the realistic neural system 
works. 
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